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Abstract 

We show that the mass-segregation solution for the steady state distrib ution of stars around a massive black 
hole (MBH) has two branches: the known weak segregation solution (iBahcall & Wolfl[T977 ). and a newly 
discovered strong segregation solution, presented here. The nature of the solution depends on the heavy-to- 
light stellar mass ratio Mh /Mi and on the unbound population number ratio Nh /Nl, through the relaxational 
coupling parameter A =4A^//M^ / [NlM-[{3+Mh /M^jj . When the heavy stars are relatively common (A^ 1), 
they scatter frequently on each other. This efficient self-coupling leads to weak mass segregation, where the 
stars form «oc r^"** mass-dependent cusps near the MBH, with indices =7/4 for the heavy stars and 
3/2 < <7/4 for the light stars (i.e. max(a// — ai) ~ 1/4). However, when the heavy stars are relatively rare 
(A <C 1 ), they scatter mostly on light stars, sink to the center by dynamical friction and settle into a much steeper 
cusp with 2 < a// < 1 1/4, while the light stars form a3/2<aL<7/4 cusp, resulting in strong segregation (i.e. 
max(a// — ai)~ 1). We show that the present-day mass function of evolved stellar populations (coeval or 
continuously star forming) with a universal initial mass function, separate into two distinct mass scales, ^ 1 Mq 
of main sequence and compact dwarfs, and ^ lOM© of stellar black holes (SBHs), and have A<0.1. We 
conclude that it is likely that many relaxed galactic nuclei are strongly segregated. We review indications of 
strong segregation in observations of the Galactic Center and in results of numeric simulations, and briefly list 
some possible implications of a very high central concentration of SBHs around a MBH. 
Subject headings: Galaxy: kinematics and dynamics — stellar dynamics — black hole physics 



1. INTRODUCTION 

Early theoretical stud ies of the dynamics and distribution 
of stars around a MBH (iPeeblesll 19721 IBahcall & Wo"afT976l 
[T97^l^ungfT980h were triggered by the discov ery of quasars 
(e.g. lMatthews & Sandagll963l:ISchmidtlll963h and the real- 
ization that many gala ctic nucl ei may contain a central mas- 
sive c ollapsed object (iLynden-B ell 1969; Wolfe & Burbidge 
Il970l) . as w ell as by the discovery of X-ray sources in globu- 
lar clusters dGiacconi et alj|1972l). which vyere then thought 



to be accretin g MBHs (e.g. IWvlledll970t iBahcaU & Woig 



\WW; see also Mille r & Hamiltonll2002l " for a recent reeval 
uation of this possibility). The main motivations for these 
studies were the prospect of detecting MBHs by the ob- 
served stellar density profile and by tidal disruption flares, 
and the po ssible role of tidal disruptions of stars in the growth 
of MBHs (Lightman & Shapiro 1 977'; C ohn & Kulsrucj 19781; 
ihapiro & Marchant 1978; Rees 1988) 

The renewed interest in this problem is driven by ob- 
servations of stars closely orbiting th e Galactic MBH 
dEisenhauer et al.l l2005l iGhez et aT] l2005h and th e accumu- 
lating data on their distribution and dynamics (Alexander 
[1999; Genzel et al. 2000, 2003; Schodel et al. 2003, 2007), 
as well as by the prospects of detecting gravitational waves 
(GW) from extreme mass ratio inspiral sources by future 
GW detectors (EMRIs : compact remnants inspiraling into 
MBHs, see review by 'Amar o-Seoane et akl 120071: iHopmanI 
I2OO6,) . EMRJ rates and properties strongly depend on the 
stellar density and the stellar dynamical processes within 
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(9(0. 01 pc) of the MBH, w here inspiraling sources originate 
(e.g. lHopman & Alexande r 2005. 2006a). 

Mass segregation occurs in dynamically relaxed systems. 
MBHs are naturally expected to lie in relaxed cores in sce- 
narios where the MBH is formed by run-away m ergers in th e 
extreme central density following core collapse (lReeslll984h . 
which occurs on timescales much longer then t he relaxation 
time, Tr (e.g. ISpitzed 119871; lOuinlanI 119961; iFreitag et all 
2006b c). Likewise, the extreme mass ratio targets of the 
planned Laser Interferometer Space Antenna^ GW detector 
(LISA) are expected to originate in relaxed nuclei, since LISA 
design is sensitive to GW from MBHs with mass < lO^M©. 
The observed correlation between the MBH mass M, and the 
typical velocity dispersion of the spheroid of the host galaxy, 
M, ( 7^, 4 < g < 5 (the Mj o relation, iFerrarese & MerrittI 
l2000t iGebhardt et~al] |2000 |) then imphes fliat su ch nuclei 
are dyn amically relaxed and very dense ( [Alexa nder 2007j; 
iMerritt et al.ll20()7h . This can be seen by assuming for sim- 
plicity j3 =4 (a higher value only reinforce these conclusions), 
and noting that the MBH radius of influence r/, ^ GM,/ oc 

1/2 

MJ encompasses a stellar mass of order M,, so that the 
number of stars there is Nh^Mt/Mi,, where Mj, is the typi- 
cal stellar mass, and the mean stellar density is Hh^Nh/rj^ °^ 



M. 



1/2 



The "nvE" rate estimate of strong gravitational col- 



lisions then implies that [r^) nha{GMi, / o^ f- oc M, ^^'^ . 
A more rigorous estimate shows that for the Galactic MBH 
(M. ~4 X lO^M,^. lEisenhauer et al.ll200l iGhez et an r2005). 
an archetype of LISA targets, 7}; ~ 0(1 Gyr) < tn (flie Hubble 

' See LISA mission website http ://lisa.nasa. gov| 
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time) and «/, (9(10^pc ^). The density in the stellar cusp 
near the MBH is orders of magnitude higher still (see below). 

Since T^ocM, ^^'^ , isolated nuclei with M, < 10^ M0 are pre- 
dicted to be relaxed. 

A single mass stellar system aroun d a MBH is expecte d 
to relax to a r^" cusp with a = 7/4 (B ahcall & Woil[l976h . 
This results from the fact that the gravitational orbital energy 
gained by the system when stars are destroyed near the MBH 
is conserved as it is shared and carried outward by the re- 
maining stars, E{r) ^E{r)N{< r)/TR oc ^-i^3-a^^a-3/2 _ 
pl2-2a ^ ^Qjjgj jBinnev & TremaindflQSTl) . When the sys- 
tem includes a spectrum of masses. Ml <M<Mh, the ap- 
proach toward equipartition by 2-body interactions decreases 
the specific kinetic energy of the high-mass stars, while that of 
the low-mass stars increases. As a result, the high-mass stars 
sink and concentrate in the center on the dynamical friction 
timescale Tdf^T^ (M) /M//, while the low-mass stars float out 
dScitzer 19 81. 

iBahcall & Wolfl (1 19771 hereafter BW77) approximated the 
mass segregation problem in the Fokker-Planck formalism, 
and solved for the steady state, angular momentum aver- 
aged stellar distribution functions (DFs) /m(£'). They found 
that near the MBH, the DFs can be approximated by power- 
laws, fM{E)ocE''" , and that the stellar current into the MBH, 
Qm{E), is very small, which leads to a specific relation be- 
tween the stellar mass and the logarithmic slope of the DF, 
Pm—M/AMh- In the Keplerian limit near the MBH, these 
DFs correspond to power-law density cusps, with Um = 3/2 + 
Pm- The WNll solution thus predicts a relatively small range 
of central concentrations, from aH=l for the heaviest stars 
in the populations, to — + 3 /2 for the lightest stars. We show 
below that the "zero-flow" assumption breaks down when the 
massive stars are too rare to scatter efficiently against each 
other, and instead sink to the center by dynamical friction 
against the light stars. As a result, the BW77 relation be- 
tween M and pM no longer holds, and the range of central 
concentrations far exceeds that of the BW77 solution. 

The Galactic Center (GC) provides to date the few avail- 
able observations that directly bear on the question of mass 
segregation around a MBH. The over-abundance of X-ray 
transients in the central ~ 1 pc of the GC was interpreted 
as evidence of a high central con centration of n eutron stars 
and stellar black holes (SBHs) ( Muno et all 12005); The cen- 
tral decrease in the surface density of the low-mass horizon- 
tal branch "red clump" giants was interpreted as evidence of 
the e vacuation of long-lived light ob jects by mass segrega- 
tion ( lLevill2006t ISchodel et aTll2007h . The dynamical upper 
limit on the distributed dark mass in the inner few mpc around 
the galactic MBH, Mdm/M. ~ fewxO.Ol (Mo uawad et alj 
lIOOllGillessenetri. 2008; G hez et alJl2008h . is still 10-100 
times larger than predicted by approximate theoretical esti- 
mates ( Morris 1993; Miralda-Escude & Gould 2000), Fokker- 
Planck calculations (Hopman & Alexander 2006b), or by the 
conservative drain limit (Alexander & Livio 2004), which 
places an upper bound on the maximal number of compact 
remnants that can avoid being thrown into the MBH by mu- 
tual 2-body scattering. These theoretical estimates are con- 
sistent with the upper bound derived from the observed lim- 
its on the diffuse X-ray in the central 0.7 pc (- 2 x lO'*) 
dPeegan & Navakshinll2007l) . 

Detailed analysis of observations of the GC dLevil I2OO6I ; 
ISchodel et al J 120071) suggests that the degree of segregation 
between the light and heavy stars is stronger than expected in 



the BW77 solution, although this interpretation of the data is 
not unique. Strong segregation in the GC i s further supported 
by analytic dHopman & Alexandeij|2006bl) and Monte-Carlo 
dFreitag et al.ll2006a[ M. Freitag. priv. comm.) mass segre- 
gation results (Fig. [U. These calculations do not take into 
account star formation and evolution, but instead assume for 
simplicity a non-evolving mass function that is based on a 
model of the GC's present day mass mass function (PMF) as 
an old, continuously star forming pop ulation with the "uni- 
versal" initial mass function (IMF) dAlexander & Sternberg! 
1 99^ I Alexander 2005, table 2.1). In these theoretical models 
at — 1.5, while 0C//>2 at r^O.l pc. A^-body simulations of 
stellar clusters with evolving stellar populations and a central 
intermediate mass black hole (IMBH) also pr ovide hints of 
strong mass segregation (Baumgardt et al.l l2004i J. Makino, 
priv. comm.). As we show below (Fig. HJ, these results reflect 
the fact that the relative fraction of SBHs in old populations 
is below a critical threshold needed for them to scatter effi- 
ciently against each other and maintain the weak segregation 
solution. Instead, they sink to the center by dynamical friction 
and settle into the strong segregation solution. 

This paper is organized as follows. The Fokker-Planck 
formulation of the mass segregation problem, the choices of 
boundary conditions and the various simplifying approxima- 
tions are described in ^ The physical meaning of the re- 
laxational coupling parameter. A, is discussed in ^ The A 
parameters of different stellar populations are explored in ^ 
A large grid of Fokker-Planck mass segregation calculations 
is presented and compared to simple analytic estimates in ^ 
The results and their implications are discussed and summa- 
rized in ^ 

2. THE ENERGY DIFFUSION FOKKER-PLANCK EQUATION 

The MBH dominates the gravitational potential within 
its radius of influence, where the dynamics are approxi- 
mately Kepleria n . Fo llowing the approach introduced by 
'Bahcall & Wolf (1976, T977) and its implementation by 
Hopman & Alexander (2006b), we calculate the approximate 
steady state DFs of a non-evolving stellar population in a 
spherical Keplerian system around a fixed MBH. This is 
done by solving the time and energy-dependent, angular 
momentum-averaged particle conservation equation that de- 
scribes the 2-body diffusion of stars in energy from a fixed 
unbound reservoir into t he MBH sink (resonant relaxation, 
iRauch & Tremain^l 19961 is not taken into account). One of 
the stellar mass groups, of mass M-^, is taken to be the refer- 
ence star It is assumed that unbound stars of mass M with 
velocity dispersion and specific energy e = \f/ ~ < 
(\ir = —(j) is the negative of the gravitational potential (j) and 
e > for bound stars), have a Maxwellian DF, 

/m(£) = Sm ^^^^2 )3/2 ^M£/<yii) , (e<0), (1) 

where n* is the number density normalization of the reference 
star. Since the number density of unbound stars 




the quantity 5^"* can be interpreted as the total density of 
stars of mass M far from the MBH, where y/^m ^ 1 ^"d 
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almost all stars are unbound to the MBH. In that case Sm is 
the asymptotic number density ratio of star M relative to the 
reference star (5* = 1 by definition). The asymptotic popula- 
tion ratios between stars of mass Mh and Ml are designated 
below by Nh/Nl = Sh/Sl- It is assumed that the asymptotic 
PMF, and its corresponding population ratios Nh/Nl, depend 
only on the IMF, on the star formation history and on stellar 
evolution physics, and are independent of mass segregation 
dynamics, which are effective only near the MBH. 

Stellar dynamics near the MBH are described here in units 
where G—1, mass is measured in units of the mass of the 
reference star, M*, specific energy in units of its velocity dis- 
persion e* = C7^, and time in units of its two-body relaxation 
time at the radius of influence. 



(3) 



where the Coulomb term is estimated as A = M,/M*. Phase 
space density is expressed in units of /* ^n^/ {Inol)^/^ and 
distance in expressed units of the MBH radius of influence 
r* = GM,/a^. In these units, the dimensionless specific or- 
bital energy is defined as x = e/ff^ =r*/2fl (a is the semi- 
major axis), the dimensionless time is defined as T = t/t-^, and 
the dimensionless DP as gM~fM/f*- 

The evolution of the dimensionless DF of each stellar mass 
group, gM, is described by the energy flow integral Qm, which 
expresses the change in energy due to 2-body scattering, and 
the angular momentum-avera ged effective loss-cone term Rm 
(iBahcall & Wolfl[T976l[T977h . 



d_ 
dx 



dx 



The flow integral is given by 



Qm{x)=Y,MM' 



M' 



■l"D 



QM{x,T)~RMix,T). 



dx' 



(4) 



-~ max(x,x') 



3/2 



. .dgM'{x') M' ,.^Sm{x) 

gM [X) — —gMi [x — 

dx M ax 



(5) 



The averaged effective loss-cone term in the diffusive 
regime (defined here a s x > x^iff =10) is approximated by 
jHopman & Alexandedl20 06b ) 



Rm[x) 



gM{x) 



Tr{x)\n[Jc{x)/Jlc] 



(6) 



where Jc{x) is the maximal (circular) angular momentum for 
energy x, 7/^ is the angular momentum of the loss-cone (taken 
here to be Jlso—^GM,/c, the angular momentum of the last 

stable orbit for | e | <C c^), and ^ \Y,m Sm (x)M^] is the ap- 
proximate 2-body relaxation time. The contribution of the 
loss-cone in the full loss-cone re gime, which typically cor- 
responds to x< 1/2, is small (e.g. Sver & Ulmel^ l 19991). and 
is therefore neglected here. Bahcall & Wolf ( 19761 119771) ar- 
gued that the loss-cone term has only a small effect on the 
energy dependence of the DF, and neglected it in their multi- 
mass models. In order to allow direct comparison of our cal- 
culations with their results, we omit here the loss-cone term 
as well, unless stated otherwise. 

Equation (HJi is integrated in time from an arbitrary initial 
DF until steady-state is achieved, subject to the boundary con- 
ditions that the DF falls to zero at some very high energy xd 



where the stars are destroyed, and that the unbound stars are 
replenished from a Maxwellian reservoir. 



gM{x>XD)=Q, gM(-x<0)=CMexp[(c7f/c7^)x] , (7) 

where the constant Cm is related to the population ratio Sm by 
(see Eq.[Tli 

Cm = (c*/ Om)^ Sm ■ (8) 

Specifically, if the different unbound mass groups are 
in equipartition (as assumed by BW77), then Cm — 
{M /Mi/f'I'^SM- Conversely, if the unbound stellar population 
has achieved equilibrium by violent relaxation (Lvnden-BelJ 
119671) . then all mass groups are expected to have similar ve- 
locity dispersions, and Cm = 5^ = 5'*. Here we assume violent 
relaxation boundary conditions^, since these are suggested 
by the lack of an observed correlation between scale height 
and sp ectral type of old stars in the Milky Way ( Kroupa et al] 
Il993h . by the absence of an observed vertical colo r gradient 
in edge-on disk galaxies dde Griis & Peletiei]|200Q) . and by 
the wea k observed spatial color gradients in elliptical galax- 
ies (e.g. iBoroson et al.lll983l:ICohenlll98 6'). We confirm the 
conclusion of BW77 that the steady state DFs do not depend 
strongly on these specific choices of boundary conditions (cf 
comparison with BW77 equipartition models in Fig. |4|. 
The spatial density profile that corresponds to the DF is 



riMir) 



gM 



(x) \J r-^jr — xAx . 



(9) 



3. THE RELAXATIONAL COUPLING PARAMETER 



As shown below in ^ the nature of the mass segregation 
solution, weak or strong, is determined by one parameter, 
which expresses the relative strength of Mh-Mh and Mh-Ml 
interactions. These can be quantified by the corresponding 
diffusion coefficients, D(Av||), D(Av|) and D(Av^), which 
enter the Fokker-Planck equation (e.g. iBinnev & Tremain^ 
Il987l Eq. 8-68). The two quadratic scattering coefficients 
are similar up to an order unity factor, and are approximated 
here as equal, D(Av|) ~ D(Av^) = D(Av^)/2. In a 2-mass 
system we represent these coefficients for brevity by the no- 
tation DsTp where S ^ L,H is the scattering star (light or 
heavy), T = L,H is the target star, and p = s,f is the process 
(scattering, D(Av^), or dynamical friction, D(Av||)). There 
are 8 possible permutations of the diffusion coefficients, 
({L,H} X {L,H} X {s,f}), however, since Dsts °= 2NsMj and 
DsTf °^ NsMs{Ms + Mt), there are only 4 distinct combina- 
tions of the mass and number ratios^. 

The relaxational coupling parameter A describes the com- 
petition between the self-coupling of the heavy stars and the 

^ Eqs. [s) generalize the definition of gM^x) and correct an er- 
ror in tlie relation between Cm and the population ratio, which ap- 
peared in Hopman & Alexander (2006b). The density profile presented in 
HoDman & Alexander (2006b) closely approximates the violent relaxation 
solution (cf Fig. [T] here), and not, as stated there en'oneously, the isother- 
mal solution. 

^ The quadratic scattering terms do not depend on the mass of the scattered 
target star because acceleration by a given gravitational force is equal for 
all masses. The quadratic terms express a random walk in velocity due to 
interactions with field stars on a timescale of ~ v-/Z)(Ai'^). In contrast, 
the first order dynamical friction term expresses the back-reaction of the field 
stars to the target star, and therefore does depend on its mass. The first-order 
term expresses a systematic deceleration on a timescale of Tdf ~v/D(Av||). 
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Fig. 1. — A Fokker-Planck model of the GC assuming a non-evolving PMF based on the continuous star-forming model shown in Fig. (2), with violent 
relaxation boundary conditions, a loss-cone term IHopman & Alexa nder 2006b) and assuming M, = 4 x 10* Mp, , = 1 M(:; , = 2 pc, and a total stellar mass 
M, within r,^. Left: The spatial number density of the 4 mass components (Eq.|9). Also shown are the best-fit logarithmic slopes aro und r = 0.1 pc. Right: Th e 
corresponding enclosed stellar number and total stellar mass as function of radius. The SBH number is consistent with the drain limit li Alexander & Livicl2004l) . 
also shown. The enclosed number of SBHs within 0.1 pc (0.01 pc) is approximately 1.2 x 10^ (1.1 x 10"), of neutron stars 3.2 x 10^ (12), of main sequence stars 
2.7 X 10" (9.4 X 10^) and of white dwarfs 2.3 x 10^ (76). 

light-heavy coupling in terms of global properties of the sys- 
tem, the mass and number ratios. 



A = 



DHHs+DHHf NhMI 



Dlh.s + DLHf NlMI 3 + Mh/Ml 



(10) 



The definition of A can in principle be generalized to a 
multi-mass stellar population by specifying the light/heavy 
mass boundary, M^^u and performing a weighted integration 

over the mass function'*. However, as shown in ^ below, 
an evolved stellar population (coeval or continuously star- 
forming) such as is expected near a MBH, is well approxi- 
mated by a 2-mass system for this purpose. 

3.1. Asymptotic mass segregation limits 

The weak segregation limit (the Bahcall-Wolf solution). In 
the limit A-^oo, which is the zero-flow (2m ^0) limit (®, 
the heavy stars dominate the population and relax to the single 
mass cusp a// = 7/4 (pn = 1/4). The light stars heat by scat- 
tering against the effectively infinite reservoir of heavy stars 
and diffuse to higher energies, thereby settling to a flatter cusp 
with ctL^ 3/2 (pL^O). 

The strong segregation limit. In the limit A ^ and when 
Mh ^ Ml, the light stars behave as a single mass popula- 
tion with aL = 1 /'^ (pl = 1 /4). The rare heavy stars sink to 
the center by dynamical friction against the effectively infi- 
nite reservoir of light stars. Their steady state distribution can 
be simply derived by noting that the dynamical friction force 
is approximately F °^ Pl/v^ °^ r^~"'-, wh ere r^^l^ is the 
Keplerian velocity (IChandrasekhari[T943h . The torque on the 
heavy star (assuming a circular Keplerian orbit, Lo^ r^/^) is 
then L^Fr^ f^-"-L oc r^^l'^r, so that r^ yili-o.^^ \^ steady 
state, the heavy star current, Nh, is independent of radius, 
N}i°^rr^n}i(r)~Q.owsX, so that the number density of the heavy 
stars must scale as n//ocr^(^/^^"^'ocr^"". It then follows that 
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E.g. A= 4(iVM2)^/ [?,{NM-)^+(Mh) (NM)i), where for a normal- 



ized PMF dW/dM, (NM'')^^ = 



/^«^M(dW/dM)dM//, 



{dN/dM)dM. 



M* (dA'/dM) dM, and {Mh) ■ 



Fig. 2. — The PMF and the relaxa tion coupling parameter for the 
Miller-Scalo IMF (Miller & Scalol [T979I) (top line) assuming 10 Gyr of 
continuous star formation (middle line) or a 10 Gyr old coeval population 
(bottom line). The IM F was evolved using the Solar metallicity stellar 
tracks of Sch aller et alj {1992) and the zero-age main sequence progenitor 
masses were mapped to compact remnants by 0.8 < M/M,;, < 1 .5 — > O.6M0 
WD. 1 .5 < M/M., < 2.5 ^ 0.7M. WD, 2.5 < M/M;, < 8.5 ^ 1 . 1 Mp WD 
8.0< M/M;, <30.0^ 1.4M. NS, 30.0<M*/M,j fOMr, SBH 1 Alexander! 
120051) . A was estimated for a light/heavy mass boundary Mi^j^ = 5M., (the 
exact choice of M^iu does not affect the value of A for the evolved PMFs). 



Equivalently, this result can be obtained by expressing the 
Fokker-Planck equations explicitl y in terms of th e diffusion 
coefficients as function of energy jLightman & Sh apiro 197"^ 
Eqs. 51a,b), and solving them under the assumption that only 
the drift term (dynamical friction) contributes. In practice, it 
is unlikely that physical systems actually have A ^0.1 (®; 
a// = 1 1 /4 should be considered an upper limit on the loga- 
rithmic slope of the number density distribution of the heavy 
stars. 

4. THE PRESENT DAY MASS FUNCTION 

The value of A (Eq. [TOb depends on the population's 
PMF. So-called universal IMFs, which extend all the way 
from the brown d warf boundary Mi ^ 0. 1 to M2 > 1 00 
(e. g. the ISalpeter.. 1955 I MF and its subsequ ent refinements, 
the lMiller& ScalollI979l and iKroupalllOOll IMFs), result in 
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Dimensionless energy x (x=10 r- 0.1 pc in GC) 



Fig. 3. — The steady state DF (on an arbitrary logarithmic scale) and 
its local logarithmic slope away from the boundaries, for M/y /Ml = 10 and 
Nh /Nl = 1 (A= 12.9). The loss-cone term is neglected. The boundary condi- 
tions are a sink at xq = 10'* and violently relaxed (equal velocity dispersion) 
unbound stars. The asymptotic BW77 solution is indicated by the dotted lines 
(Pi. = l/4().p« = 7/4). 

evolved populations (old coeval star-bursts or continuously 
star forming populations) that naturally separate into two 
mass scales, the ^ IM© scale of low-mass main-sequence 
dwarfs, white dwarfs and neutron stars, and the ~ lOMp scale 
of stellar black holes, and typically have A < 0.1 (Fig.|2|i. Such 
evolved populations are well-approximated by the simple 2- 
mass population model. 

Generally, 10 Gyr old, continuously star-forming popu- 
lations with a power-law IMF, d/V/dMocM"^, have A< 1 
for 7> 1.8, and A> 1 for 7 < 1.8 (assuming Ml = Q.IMq, 
Mh = I2OM0, Mi^fH = 5Mq). This impHes that even flatter 
IMFs than the Salpeter IMF (7=2.35) lead to strong mass 
segregation. 

It is not clear what is the typical PMF, and the corre- 
sponding value of A, in galactic nuclei. There are in- 
dications that star formation deep in the potential well 
of a MBH can be v e ry dif f erent from that in the field 
jLevin & Beloborodov '2003': jMilosavljevic & Loebl l2004t 
iPaumard et al. 2006; Levin For example, the volume- 

averaged stellar population in the central few pc of the 
GC is reasonably well approximated by a 10 Gyr old, 
continuously star-forming pop ulation with a universal IMF 
([Alexander & Sternberdl 19991) and has A < 0.1. In contrast, 
a recent analysis of late type giants in the inner 1 pc of 
the GC suggest s that the IMF there could be a flat 7^0.85 
power-law (Maness et alJ i2007). This corresponds to A > 10 
in the inner ~ 1 pc, if that region evolved as an isolated sys- 
tem. However, that is unlikely, since the PMF reflects not only 
the local star formation, but also the volume-averaged popula- 
tion in the inner few pc, which contributes unbound stars and 
stars that diffuse in phase space from lower to higher energies. 
It is therefore quite possible that in some galactic nuclei the 
weak mass segregation solution is realized, while in others the 
strong mass segregation, depending on the system-to-system 
scatter in the PMF. 

5. RESULTS 

Figure (O shows an example of a steady state solution of the 
DFs calculated for a 2-mass system in the weak segregation 
regime with Ml/Mh ~ 10 and Nh/Nl = 1 (A~ 13), neglect- 
ing the loss-cone term. While the logarithmic slopes of the 
DFs are not constant, they vary only slightly away from the 



boundaries at x = and xo- 

When comparing the DF with observations of stars in the 
GC, a relevant energy scale is x~ 10, which corresponds to 
orbits with a semi-major axis of a = r*/2x = 0.1 pc for r^, = 2 
pc in the GC. On that scale the cusp is Keplerian, but there 
are still enough obs erved stars for meaningful statistics (e.g. 
ISchodel et al.ll2007h . Figure (|4]L) shows the local logarith- 
mic slopes of the DFs, pi and pn, at x= 10, as function of A, 
for the same mass ratios modeled by BW77, Mh/Ml = 1.5, 
3 and 10. The transition between the weak and strong mass 
segregation solutions at A~ 1 is clearly seen. It is interesting 
to note that the set of models calculated by BW77 all happen 
to lie at A> 1, which explains why they found that ". . . there 
are no dramatic changes in the shapes of the curves, despite 
the fact that the unbound population ratio was varied by a fac- 
tor of 64.", and wrongly concluded that the weak segregation 
solution is a universal solution. 

A useful average of the logarithmic slope is obtained from 
the corresponding stellar density curve (Eq. |9l cf Fig. [T]), 
which can be directly related to the observable stellar density 
distribution. A comparison of Figs. (|4]L) and (|4]R) shows 
that the shape of the curves is hardly affected by the method 
chosen for deriving the logarithmic slopes. We also repeated 
these calculations with a loss-cone term (Eq. [H, and con- 
firmed the BW77 conclusion that the effect on the shape of 
the DF curve is small, ther eby j ustifying its omission. 

The asymptotic limits ( ^3.11 ) are also clearly seen in Fig. 
(01). When A^oo, the weak segregation BW77 solution holds, 
with ph = 0.25 and pi ^ Ml/4Mh- When A^O, the light 
stars, which dominate the population, assume the single mass 
population DF, pi =1/4, while the heavy stars concentrate to 
the center. For low mass ratios, Mh /Ml < 4, where dynamical 
friction is less efficient, the heavy stars approximately obey 
the BW77 relation, p// = {Mh /Ml) pl ^Mh /'\-Ml. For higher 
mass ratios, Mh /Mi>4, the heavy stars approach the dynam- 
ical friction limit, pn ^5/4. 

The transition between the weak and strong mass segrega- 
tion solutions is a reflection of the breakdown of the zero-flow 
assumption as A ^ 0. The dimensional scale /* of the dimen- 
sionless stellar current into the MBH, Qm{£) =lM{£-)/h, is 
/* ^ (9(A'm* (< (note that when the loss-cone is ne- 
glected, the steady-state current is independent of energy). 
This is also the order of magnitude of the stellar current into 
the MBH in a single mass po pulation out of equilibrium. 
However. iBahcaU & Woi3 ([197 6. Eq. 63) show that, neglect- 
ing the loss-cone, the steady-state current in a single mass 
population is Qm ~ 8 /xd ^ 0, and that steady state zero-flows 
further ir nply that ph(£)/Mh — Pl{£) /Ml in multi-mass pop- 
ulations ( lBahcaU & Woigil977L Eqs. 41-46). The relevant 
physical timescale for the heavy star current 2// is not the Ml- 
Ml interaction timescale 1/f,^ = F/, oc Dlls + Dll/ ~ ANlM^ 
(here we choose ^Ml), but rather the rate of Mh-Mh and 
Ml-Mh interactions, Fh °= Dhhs + Dlhs + Dnnf + DLHf ~ 
NlMI [4 + Mh/Ml + ANhM^j/NlMI] . To compare meaning- 
fully the heavy star current Qh across the range of our mod- 
els, we rescale it to 1 /Fh instead of f*, and normalize by the 
number of heavy stars to obtain the current per star, 

^" ^"Nh Fh ^" Nh / V Ml ^ NiMI ) ' 

(12) 

Figure (|5]) shows Qh as function of A for different mass ra- 
tios. For A> 1, which is the range explored by BW77, we 
confirm their zero-flow result. In contrast, for A<C 1, we find 
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Fig. 4. — The power-law indices pi n as functions of A for various mass ratios, assuming violent relaxation boundary conditions, and neglecti ng the loss-cone 
term. Left: The indices, as derived from the logarithmic slopes of gi n at x= 10, compared to the logarithmic slopes of the DFs calculated bvlBahcall & Wolll 
09771) . Right: The indices, as derived from the average logarithmic slopes of the densities 3/2 + pt n in the ran ge r = 0.003-0.3 pc (for r^ = 2 pc) with the 
expected pn for a model of the GC assuming conti nuous star formation with the Miller-Scalo IMF (Fig. l2l iAlexander & Stemberelfl999l : IAlexandeni2005i) , and 
for a mass-function model of globular cluster M15 IMurphv et al.ll997h . assuming it harbors an IMBH. 

Figure (|4]R) shows the expected segregation of the massive 
stellar objects in the GC (based on the old, continuous star- 
formation PMF model. Fig. |2]i, and in globular cluster M15, 
assuming it harbors an IMBH (based on the tentative PMF 
model of Murphv et al. 1997, which assumes that the most 
massive remnants in the cluster are 1 AMq, neutron stars). The 
GC is expected to lie in the strong relaxation regime (A = 
0.05), while M15 in the weak relaxation regime (A^ 18). 

We refrain here from a more detailed quantitative analysis 
of the numerical results because the exact values of pl h{^) 
depend somewhat on the way these are evaluated (locally, or 
from the density curve), and because the convergence of our 
numerical Fokker-Planck solver becomes progressively worse 
in the limits A— ^ and Mh/Ml ^ 1. However, any numeri- 
cal inaccuracies in those limits are unlikely to have practical 
implications, since realistic stellar systems are not expected to 
have A<0.01 (a possible exception could be the hypothesized 
steady flow of IMBHs from dense stellar clusters to the MBH, 
iPortegies Zwart et alj|2006l) . 
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Fig. 5. — The normalized cun'ent of the heavy stars, Qh (Eq. I12t , as 
function of the relaxational coupling parameter, A. The horizontal line marks 
Qo = ^Ixd, the zero-flow in a single mass population (Bahcall & Wolf 197(3). 
For AS>1, Q/f~Qo~0. However, for A ^ 1 in the strong mass-segregation 
regime, Qh 0( 1 ), indicating that massive stars flow inwards at the maximal 
possible rate, the dynamical friction rate. The transition between the two flow 
regimes is also apparent in the energy-averaged currents in models with a 
loss-cone term. 

that the zero-flow assumption no longer holds, Qh ^0(1), 
and the heavy stars sink to the MBH at the maximal possi- 
ble rate, the dynamical friction rate. Figure (|5]l also shows 
that the physically required inclusion of a loss-cone term in 
the Fokker-Planck equations (Eq. |6]l increases the current, but 
the energy averaged flow in the weak segregation limit A^ 1, 
while not quite as small as without a loss-cone, is still strongly 
suppressed^. It then follows that the relation ph/Mh^Pl/Mi 
is approximately satisfied in this limit even with a loss-cone, 
as found by BW77 and confirmed here. In the limit strong 
segregation limit A^ 1, the dynamical friction-driven current 
of the heavy stars is insensitive to the presence or absence of 
a loss-cone, since T^i^ {Mi/ Mh)Tr. 

^ In the presence of a loss-cone, the cun'ent in the weak segregation 
limit is set by the diffusion rate into the loss-cone at xjjff (©, Ih ~ 
A'«(-Vdiff)/['Hlog([7c(jCdiff)Mc-)], where /a~r^' is the (M;^ -dominated) re- 
laxation time (cf Eq. |6). It then follows that Qh (^diff ) ~ ("'diff) / Q^h^h ) ~ 
[1 /(2W]^"'"' / log 0-* )/ ^32^ ~ 0.01 for .Vdiff = 10. 



6. DISCUSSION AND SUMMARY 
6.1. Strong segregation and other relaxation processes 

Strong mass segregation, or mass segregation instability, 
shares some common featur es with the S pitzer, or equiparti- 
tion, instability in a cluster ('S pitzeiiri969 ). where the heavy 
stars decouple from the light ones and evolve away from 
equipartition by forming a dense sub-system in the center with 
a much higher velocity dispersion, <7^/of >Ml/Mh- How- 
ever, these two instabilities are distinct effects. The Spitzer in- 
stability occurs only when the heavy stars are relatively com- 
mon in the population, NhM^J^ /NlMI^^ > j3, where j3 ~ 0. 16 
(for Mh :>Ml and NrMh -^NlMl). In conti-ast, in the Keple- 
rian potential near a MBH, the Jeans equation dictates that 
alJal = (5/2 + pl)/ {5/2 + pn) ^ I (Alexander & KumaJ 
I2OOII) . and so equipartition is never achieved, irrespective of 
the heavy-to-light mass or number ratios. Strong mass segre- 
gation is an instability in the spatial distribution of the heavy 
stars, which occurs in the opposite limit to the Spitzer instabil- 
ity, when the heavy stars are relatively rare in the population, 
A ryNHMH/NiMnK l (for Mh :$>Ml, Eq. [TO]). 

iPerets et akl (l2007h focused on the relaxation of light ob- 
jects, stars, by heavy objects, massive perturbers (e.g. giant 
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molecular clouds or clusters, with Mh ^ Mi), and expressed 
the efficiency of massive perturber-induced relaxation relative 
to star-star relaxation by the parameter ji2 =NHMjj/NLMj^ 

{DHLs +DHLf)/ {DLLs +DLLf) ^ 3NhM]j / ANlMI) . When 
/X2^ 1, massive perturbers dominate stellar relaxation. Here 
we focus on the dynamics of the heavy objects, the SBHs, 
and so the relaxational coupling parameter A is defined to also 
take into account the interactions between the heavy masses. 
The parameter 112 addresses the question "Which mass com- 
ponent dominates the relaxation of the light stars?", while A 
addresses the question "Which mass component dominates 
the dynamics of the heavy stars?". The two are related by 
A = I124/{3+Mh/Ml). 

Our approximate treatment of the mass segregation process 
neglects the relaxation of angular momentum to near radial 
("loss-cone") orbits. This, rather than diffusion in energy, is 
the primary channel for stellar destruction by the MBH. A 
full treatm ent of the m ass segregation problem in {E,J) phase 
space (e.g. Cohn & K ulsrud 1978) is beyond the scope of this 
work. Here we follow BW77, who treated the problem ap- 
proximately in E only, and who further showed that the ne- 
glect of an effective loss-cone term (Eq. |4]i did not much 
change the shape of the DFs. We confirmed that this con- 
clusion also holds for our mass segregation models (Fig. [Sj. 

6.2. Possible implications of strong mass segregation 

Strong mass segregation occurs in stellar systems with a 
relatively lower fraction of SBHs, that reach a higher central 
concentration of SBHs very close to the MBH, compared to 
systems with a higher fraction of SBHs that undergo weak 
segregation. To compare two such systems, and to determine 
which has more SBHs enclosed inside some given volume 
around the MBH, it is necessary to specify the comparison 
procedure (e.g. assuming an equal total stellar mass; or the 
same MBH mass and the M,/a relation; or an equal num- 
ber of SBHs within r/,). The choice depends on the question 
of interest. Here we do not address such quantitative issues, 
but limit ourselves to briefly listing some processes that are 
expected to be affected by the degree of segregation. 

Accelerated relaxation. The degree of mass segre- 
gation affects both the non-coherent 2-body relaxation 
timescale, which scales as 1 //M^(dA^/dM)dM (see ®, 
and the resonant relaxation timescale, whic h scales as 
asJ_M(dA^/dM)dM / j M"^ {dN / dM)m (Rauch & Tremaind 
Il996l) . In particular, the stronger the mass segregation, the 
shorter is the resonant relaxation timescale, which does not 
depend on the number of stars, but only on their typical mass. 
Efficient resonant relaxation near the MBH may affect stellar 
orbits and accretion disk dynamics there (see below). 

GW event rates. The GW EMRI rate is determined by the 
number of potential GW sources within the critical radius 
rcrit, which demarcates the boundary between compact object 
that inspiral into the MBH those that plunge (infall) into it. 
The critical radius is a function of the relaxation time, and 
to good approximation the EMRI rate is T ^ N[rcrit{TR)\/TR 
dHopman & Alexandeij|2005,) . Strong segregation will affect 
the EMRI rate both by modifying the 2-body relaxation time 



and by affecting the number of stars enclosed inside rcrit, as 
well as by decreasing the resonant relaxation timescale. Sim- 
ilarly, the rates of detectable GW bursts from fly-bys near the 
Galactic MBH strong ly depend on the number of SBHs near it 
fflop man et alj2007l) . Strong segregation may also affect GW 
emiss ion from close SBH -SBH interactions in a very dense 
cusp dO'Learv et alj|2008l) . 

SBH-star interactions. A higher central concentration 
of SBHs affects the probability of SBH-star interactions, 
which can lead to the ran domization of stel lar orbits, the 
heating of a stellar disk dPere ts et al.! |2008l) , the 3-body 
exchange capture of ma ssive young stars near the MBH 
( [Alexan der & Livid 120041). or the ejection of hyper-velocity 
stars dP'Leary & Loeb 2008). 

SBH-accretion disk interactions. A higher central concen- 
tration of SBHs within lO'^^ gravitational radii of the MBH 
could exert coherent torques on the accretion disk, warp it and 
possibly affect its hydrodynamics (Bregman & Alexander, 
2008, in prep.). The SBHs may shock the disk as they cross 
it, experience drag by it, and be carried by it to the MBH (e.g. 
Artvmowicz et al.J993;.Navakshin et ali2004 ; Subr & KarasI 
20051k 

Enhanced gravitational lensing. SBHs projected near the 
Einstein angle of the MBH can strongly modify the gravi- 
tational lensing properties of the MBH, in a way similar to 
the effect of a plane t orbiting a Galactic star that is lensing a 
backg round source ( [Alexander & Loebll2001l; [Chaname et alj 
l200l . 

6.3. Summary 

We show that the steady state solution of a relaxed multi- 
mass stellar system around a MBH has two branches: the 
known weak (Bahcall-Wolf) mass segregation solution, where 
the difference in the degree of central concentration of the 
light and heavy stars is relatively small, and a newly discov- 
ered strong segregation solution, where the difference is much 
larger. The nature of the solution is determined by the global 
properties of the system (the mass ratio between the heavy 
and light stars, Mh/Mi, and their number ratio far from the 
MBH, Nh/Nl) through the relaxational coupling parameter, 
A^ANnMjj / [NlMI{3+Mh / Ml)] . Strong mass segregation 
occurs when the heavy stars are relatively rare in the popu- 
lation (A<C 1), and sink to the center by dynamical friction. 
Weak mass segregation occurs when the heavy stars are com- 
mon in the population (A^ 1) and settle to the single mass 
stellar cusp solution. We show that relaxed old coeval or con- 
tinuously star-forming populations with a universal IMF typi- 
cally have A<0.1, and thus settle to the strong mass segrega- 
tion solution around a MBH. 
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